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TECHNICAL MEMORANDUM 


A “KANE’S DYNAMICS” MODEL 
FOR THE ACTIVE RACK ISOLATION SYSTEM 

1. INTRODUCTION 


The vibratory acceleration levels currently achievable, without isolation on manned space 
structures, exceed those required by many space science experiments. 1 ^ Various active isolation devices 
have been built to address this need. The first in space was called Suppression of Transient Accelerations 
by Levitation (STABLE), which uses six independently controlled Lorentz actuators to levitate and 
isolate at the experiment or subexperiment level. 5 It was successfully flight tested on STS-73 (United 
States Microgravity Laboratory-02 (USML-02)) in October 1995. Building on the technology 
developed for STABLE, Marshall Space Llight Center is developing a second-generation experiment- 
level isolation system: GLovebox Integrated Microgravity Isolation Technology (g-LIMIT). 6 This 
compact system will isolate microgravity payloads in the Microgravity Science Glovebox. 

A second experiment- level isolation system, the Microgravity Isolation Mount (MIM), was 
launched in the Priroda laboratory module, which docked with Mir in April 1996. 7 MIM uses eight 
Lorentz actuators with centralized control. It has supported several materials science experiments since 
its implementation in May 1996. A modified version of MIM (MIM II) supported additional experiments 
on STS-85 in August 1997. 

Boeing’s Active Rack Isolation System (ARIS), in contrast to the above payload isolation sys- 
tems, has been designed to isolate at the rack level. An entire international standard payload rack (ISPR) 
will be isolated by each copy of ARIS on the International Space Station ( ISS ). The risk mitigation 
experiment for ARIS was conducted in September 1996 on STS-79. 8 Each of ARIS’s eight electrome- 
chanical actuators requires a two rigid-body model. When the ISPR (“flotor”) is included, the total 
isolation system model contains 17 rigid bodies. 

In order to provide effective model-based isolation, the task of controller design requires prior 
development of an adequate dynamic (i.e., mathematical), model of the isolation system. This technical 
memorandum presents a dynamic model of ARIS in a state-space framework intended to facilitate the 
design of an optimal controller. The chosen approach is the method of Thomas R. Kane (Kane’s 
method). 9 The result is a state-space, analytical (algebraic) set of linearized equations of motion for 
ARIS. Interim versions of this work appeared as references 10 and 11. 



2. THE CHOICE OF KANE’S METHOD 


There are fundamentally two avenues for deriving system dynamical equations of motion: vector 
methods and energy methods. Both avenues lead to scalar equations, but they have different starting 
points. Vector methods begin with vector equations proceeding from Newton’s laws of motion and 
energy methods begin with scalar energy expressions. The former category uses approaches built around 
(1) momentum principles, (2) D’Alembert’s Principle, or (3) Kane’s Method; and the latter uses 
approaches built around (1) Hamilton’s Canonical Equations, (2) the Boltzmann-Hamel Equations, 

(3) the Gibbs Equations, or (4) Lagrange’s Equations. 

Although some problems might lend themselves better to solution by other approaches, Kane’s 
method appeal’s in general to be distinctly advantageous for complex problems. As a rule, of the above 
approaches, those that lead to the simplest and most intuitive dynamical equations are the Gibbs Equa- 
tions and Kane’s Equations. Of those two approaches, the latter is the more systematic and requires less 
labor. The reduction of labor is particularly evident when one seeks linearized equations of motion, as 
proved to be necessary in the present case (due to the otherwise excessive algebraic burden). 

An over-view of Kane’s approach to developing linearized equations of motion is presented in 
reference 10, along with a summary of the relative advantages of the method. See Kane and Levinson 
for more extended treatments. 9 ’ 12 
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3. DESCRIPTION OF ARIS 


The total dynamical system S consists of the stator S ( ISS and the integral frame, from the mo- 
tion of which ARIS isolates the ISPR), the flotor F (the ISPR), eight electromechanical actuator assem- 
blies, and the umbilicals (see fig. 1). The flotor is connected to the stator by the eight actuator 
assemblies and by a variable number of umbilicals. The actuator assemblies also (and fundamentally) 
act as the vibration isolation devices. 
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Each actuator assembly consists of a Lorentz (voice-coil) actuator, an arm, an upper stinger, a 
push-rod, a lower stinger, and a position sensor. (See fig. 2 for a kinematic diagram and fig. 3 for a 
computer-aided design (CAD) drawing of a single actuator.) One end of each actuator arm is connected 
to the flotor through a cross-flexure, which allows the flotor a single rotational degree of freedom (DOF) 
with respect to the stator. The other end of the arm is connected to one end of the push-rod through the 
upper stinger, a wire of very high torsional stiffness. Each upper stinger provides two rotational DOF’s 
in bending. The opposite end of the push-rod is connected to the stator through the lower stinger, another 
short wire that allows three rotational DOF’s (two in bending, one in torsion) with respect to S. 


Lorentz Coil 



Figure 2. Kinematic diagram, including the /th actuator assembly and the umbilical. 
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Field Assembly 


Coil Assembly 



Position — 1 
Detector 


3-DOF Hinge 

and Launch Sleeve 


Figure 3. A single ARIS actuator. 


Each stinger is modeled as a massless spring. The umbilicals are also considered to be massless; 
they are modeled together as a single, parallel spring-and-damper arrangement, attached at opposite ends 
to stator and flotor, at effective umbilical attachment points S u and F u , respectively. This effective 
umbilical applies both a force and a moment to the flotor. The force is assumed to act at point F u . 

The stator, the flotor, and each actuator arm and push-rod are considered to be rigid bodies with 
mass centers at points S*, F , A*, and P* , respectively. The superscript * indicates the mass center of the 
indicated rigid body; the subscript i corresponds to the /th actuator (/ = 1, ..., 8). All springs (cross 
flexures and stingers) are assumed to be relaxed when the ISPR is centered in its rattlespace (the “home 
position”). 
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4. COORDINATE SYSTEMS 


With the ISPR in the home position, fix eight right-handed, orthogonal coordinate systems in the 

flotor, one at each of the cross-flexure centers. Let the /th coordinate system (/ = 1, ..., 8 ) have origin F ; 

located at the center of the /th cross flexure, with axis directions determined by an orthonormal set of 

unit vectors f \ j = 1, 2, 3). (The overhat indicates unit length, the index / corresponds to the /th actuator 
J '•i 

assembly, and the index j distinguishes the three vectors.) Orient the unit vectors so that /. is along the 

/th arm, toward the /th voice coil; / is directed parallel to the other segment of the /th arm and toward 

^ j *i *i 

the upper stinger (which is located at A 2 )\ and f ^ is in the direction / 1 x f 0 (along the intersection of 
the two cross-pieces of the /th cross-flexure). 

Fix a similar right-handed coordinate system a 1 : (j = 1, 2, 3) in the arm of each actuator. Locate 

J ' Aj 

each system a • such that it is coincident with the corresponding flotor-fixed coordinate system / . when 
J J 

the flotor is in the home position. 

At the respective lower stingers (points 5 ( ), place eight push-rod-fixed coordinate systems p\ 
and eight stator- fixed coordinate systems s l j. Orient these 24 coordinate systems so that when the 
stingers are relaxed (i.e., with the ISPR in the home position), the coordinate directions p l . and s) are 
coaligned for the /th actuator with p\ (along with s 2 , in the home position) directed from toward Ah. 

Finally, define a primary, central, flotor-fixed, reference coordinate system with coordinate 

/V 

directions f .. All other flotor-fixed coordinate systems f . are assumed capable of being referenced 
(e.g., by known direction cosine angles) to this system (see eq. (4)). 


5. ROTATION MATRICES 


Let the q'j coordinate system rotate, relative to the / . coordinate system, through positive 

angle q[ about the f ^ axis. Similarly, let the orientation of the a'j coordinate system, relative to the p l 

coordinate system, be described by consecutive positive rotations q 2 (about the pj axis) and q l 2 (about 

the moved three-axis). Let the orientation of the p l . coordinate system, relative to the s'- coordinate 

J • , • J 

system, be described by consecutive positive rotations q\ (about the I 3 axis), q l $ (about the moved two- 
axis), and q g (about the moved one-axis). 
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Let c'j and s l j represent the cosines and sines of the respective angles q l j . Then the rotation 
matrices among the several coordinate systems for the /th actuator assembly are as follows: 


J J 
C 4 C 5 

J J 
S 4 C 5 

1 

Co 

L/V-*- 

1 

£1 

J J . J J J 
-s 4 c 6 +c 4 s 5 s 6 

J J 1 J J J 
C 4 C 6 + S 4 S 5 S 6 

J J 
C 5 S 6 

p 

^2 

J J 1 J J J 
S 4 S 6 + C 4 s 5 c 6 

J J l J J J 
-C 4 S 6 + S 4 S 5 C 6 

J J 
C 5 C 6 

p 

1-3 J 


and 


/v / 

- 1 


4 

44 

44 

/V / 
—2 

► = 

-4 

44 

44 



0 

-4 

4 


Q.\ 



4 

0 

f 

Lx 

q 2 

► = 

-4 

4 

0 

f 

2-2 

-3 J 


0 

0 

1 

f 

2-3 


(2) 


( 3 ) 


Finally, define a rotation matrix between the eight flotor-fixed coordinate systems f . and the single 

a , J 

flotor-fixed reference coordinate system / . : 


l 


/A fk th 

A 

Ll 

f-2 

, = 

fix fh fh 

1-2 

i 3 


fix fh fh 

^3 


( 4 ) 


6. GENERALIZED COORDINATES FOR S 


The 48 angles q'j are the generalized coordinates of the system. For the /th actuator, the six 
associated generalized coordinates are as follows: q[ is the angle at the cross flexure of the /th actuator; 
q*2 and <73 are the angles at the upper stinger; and <74, <75, and q^ are the angles at the lower stinger. 
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7. GENERALIZED SPEEDS FOR S 


Define generalized speeds u l j for the system as the time rate of change of the generalized coordi- 
nates of S in the inertial reference frame: 


u) = q) (for j = 1, .... 6; i = l , .... 8) . (5) 

8. ANGULAR VELOCITIES OF REFERENCE FRAMES AND RIGID BODIES 


Designate the reference frames corresponding to the stator, the /th push-rod, the zth arm, and the 
flotor, by the symbols S , /), Aj, and F , respectively. Let Sj and F; represent, respectively, the 


coordinate systems in S and F defined respectively by s\ i 3 


and 


h Li £3 


"I T 


Two 


intermediate reference frames were introduced previously to permit describing the angular velocity of 
each push-rod relative to the stator; designate those intermediate frames corresponding to the zth actuator 
assembly by Rj and Q v Another intermediate reference frame was previously introduced between 
frames P t and A,- ; designate this by Tj . 


Let each intermediate reference frame have a frame-fixed, dextral set of unit vectors. Indicate the 
unit vectors for each of these frame-fixed coordinate systems by using the corresponding lowercase 
letter ( r‘j corresponding to R t , etc.). The following, then, give the expressions for the angular velocities 
of the various reference frames and rigid bodies of 5 : 


F i 0) A > 

in 

II 

(6) 

p ico T i 

i 

— u 2P_i > 

(7) 

r i a> A ‘ 

CO 

<■^.1 

•-co 

S 

II 

(8) 

s i oj R i 

i *i 

= U4S3 , 

(9) 

R i g)Qi 

i 

= u 5L2 . 

(10) 
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( 11 ) 


Q i(0 p i =u l 6 q[ 


Using the addition theorem for angular velocities, the angular velocities of the rigid bodies 
of S are 


and 


S; 


CO A ‘ 


i yv/ i i a. i i a,/ i a. I 

M2p| + + M4S3 + + u^q l { 


( 12 ) 


■Si- 


ft) 


P, /a./ /a,/ /a,/ 

' = M4S3 + M 5 r 2 + M 6 ^/] 


(13) 


ft)^ 


7 a./ 7 a. 7 7 a. 7 7 a. 7 

= U 2 P^ "1" ^3-3 "1" "I" ^5—2 ^6 d ] 

9. BASIC ASSUMPTIONS 



(14) 


In the subsequent development of the ARIS equations of motion, it is assumed that ARIS works 
as intended; i.e., that the ARIS controller prevents the ISPR from exceeding its rattlespace constraints. 
It is also assumed that the small-angle approximations hold for angles q). Angular velocities and angu- 
lar accelerations are assumed to be small as well. This means that the use of first-order linear perturba- 
tions will permit the full nonlinear equations of motion to be approximated accurately by a set of 
first-order linear differential equations. Finally, it is assumed that the angular velocity of the stator is 
negligible and that the stator translational velocities and accelerations are small. 


10. LINEARIZED VELOCITIES OF THE CENTERS 
OF MASS FOR THE RIGID BODIES OF S 


Represent by r AB the position vector from arbitrary point A to arbitrary point B. Define the 
following position vectors using the indicated scalars: 


r _ FiA 2 =/{#> + 4a i , 

(15) 

SjA 2 A 

r l = hP 2 , 

(16) 

S[ Pi i '•i 

L =P2P 1 . 

(17) 


9 



and 


( 18 ) 


AAi 


l A. I 7 /V 7 

a{a[ + a2a2 


JiF 


/l'/' 1 +/2?2+/3/3 • 


(19) 


First-time derivatives of the appropriate position vectors, under the stated assumptions, yield 
expressions for the velocities of the centers of mass for the 17 rigid bodies. The following expressions 
are the linearized velocities for those centers of mass. (The presubscript indicates that the expressions 
are linearized; the pre superscript indicates the reference frame assumed fixed for purposes of the 
differentiations.) 


S ‘lY P ‘ = p l 2 (- u l 4 s[ +« 6 ^) 0 = !’ •••’ 8 ) 


(20) 




/- 


-a2«3 — {@2, A (3 p j + 1 < 3 1 W 3 + r/ j ui^ j p + CI2M2. — a \ u 5 A y 1 ! ^3)^6 


( 21 ) 


and 


Si F* 
l v - = 


* 1 

Ei 


/2«l - V2«3 "(V2 + /J)«4 + v 3“5 

h “//“l ~ v 3 ( M 2 + «6 ) + V 1 ( u 3 + u 4 ) M + V 2 M 2 “ W «5 + ( v 2 + 4 ) u 6 


-1 

^3 


where 

v{=//-4 , 

and 

vWi • 


( 22 ) 


(23) 

(24) 

(25) 


11. LINEARIZED ACCELERATIONS OF THE CENTERS OF MASS 
FOR THE RIGID BODIES OF S 


Taking the time derivatives of the respective linearized velocity vectors yields expressions 
for the linearized accelerations of the centers of mass for each rigid body. Note that the linearized 
velocity vectors may be used in this step — the full nonlinear accelerations need not be determined. 
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This is a tremendous savings of effort, which would not be afforded if Newton’s Second Law were 
applied directly instead of Kane’s approach: 


S: p* i ( ■ i ■ i \ 
l a ‘ = P2[-U4S 1 +u 6 s 3 ) , 


(26) 


S \a A i = 


_ 44 — (4 T 4 )4 j4 T (44 L 44)4 


and 


+ 


44 — 44 + (4 T (3 ) 4 


-3 ’ 


V* = 


/ 2«1 - 44 - (4 + 4)4 + 44 


* 1 
£l 


(27) 


+ 


-/lM - 4(4 +4)+ 4(4 + 4) 


+ 


44 -44 +(4+4)4 


-1 

£3 


(28) 


12. LINEARIZED PARTIAL VELOCITY VECTORS FOR THE POINTS OF S 
AT WHICH THE CONTACT/DISTANCE FORCES ARE ASSUMED TO ACT 


The partial velocities and partial angular velocities are formed by inspection of the relevant 
velocity vectors. These partial velocities are then (and the order here is c racial ) linearized by neglecting 
higher order terms. 


12.1 Linearized Partial Velocities of P/ 

For the /th push-rod, the linearized partial velocities are 

* 

S jv? =0 (for r=l,2,3) , 

jfc 

S: P: i l i \ 

lU =~P2\S\ + c M^2j . 

S; P; i i ~i 

=P2%h . 
and 


(29) 

(30) 

(31) 
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12.2 Linearized Partial Velocities of A,- 


For the ixh ann, the linearized partial velocities are 


S: A: 

l l v - 1 * =0 ’ 


* 

S )Y .2 = “ 4 ^ 2 £2 + ( a 2 + ^(^3)^3 » 

S / Y 3 = ~( a 2 + a i^3 )pi + («{ - a m)p' 2 + «k2^3 . 

‘V 4 ' = ~{ a 2 + a F/3 + 4 )p 1 + («{ - a l c l l ?,)p2 ~ [ a 2 c l5 + a l<?6 + ^ 5)^3 

‘Vs' = ( a 2<?2 + a 2<?6 + 1 3<16)p\ ~ a WeP! 2 + 2^3 “ a l j/^ > 
V(i' = -ai&K + [ a 2 + a W?> +4)^3 • 


12.3 Linearized Partial Velocities of F 


For the flotor, the linearized partial velocities are 


Si F* A A 2 1 

A 1 =f 2 f l -flf 2 . 


lS| /F2 - - (^2^2 + v \)p 2 + ( V l<73 + v 2~ v Al2 ~ > 


( _v lM _ v 2 + Zl'dl )p) + (W “ v 2 c l\ + fhll )p\ + V/ 2 P 3 » (41) 

lS /ll4 = ~ v \ c l\ ~ v 2 + V \{ c l2 + %) + f\ ( l\ ~ l \ t\ 

+ (v) - v\q\ + v\q\ + f 2 q\ )p\ - [v\q\ + v\q\ + l\q\ ] Jp\ > , (42) 
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Si V F 
1—5 


and 


v\{<12 + 96 ) + v 3 + 4% P\ ~ v\<llp 2 - (W - v\q\ + fW\)l}\ . 
‘ S /^6 = _ ( v 2<?2 + v 3 )/4 + ( V M + v 2 ~ v \ c l2 ~ f\<k + 4)p 3 • 


( 43 ) 

( 44 ) 


12.4 Linearized Partial Velocities of F„ 


F F 

Define measure numbers for r “ as follows: 


F F a 1 a 1 

r u =X F f +Y F f 0 +Z F f . 

- r u — 1 r u —2 r u —3 


Since 


1 S 1 & /?* 

/V= /— r + 


du r 


Vx r F F u 


( 45 ) 


( 46 ) 


the linearized partial velocities for the umbilical attachment point F u can be expressed as follows: 

S )Yi u = S )vf +Y Fu f\-X Fu f \ , 


and 


Si F u S, f* , 
W= /-3 + 


X 


Fu 


(<0 -(4) 


l Fu 


Ei + 


X Fu +Y Fu\ 


Si F S t f* 


/-4 


X Fu (<7l “ (4) + Z Fu |<Z2 +%)“ Y Fu 


+ 


K + - ft) + - [X Pu q 6 + Y Fu q\]p l 3 , 


Si F u Sj /r , 
Z-5 = Z-5 + 


Z Fu +Y Fu\ 


(<72 + <76 ) p\ ~ [ X Fu%)p 1 2 ~ Y Fu[q 1 “</]) + X Fu 


'*> 1 
P 3 


S ) y.(“ ~ S ) 1:6 _ ( Z /'« + Y Fu c l2 ) £2 + + Z Z-« (<?3 “ <-l\ ) “ Z Fttt2 


713 


( 47 ) 


^“=^2* -(Zf«+ W)£J+[i> w +Xf«(«5-91 1 )- z F»92]pJ . (48) 


{q\-q\)Y 2 + x Fuq\p^ , ( 49 ) 


( 50 ) 

( 51 ) 

( 52 ) 
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12.5 Linearized Partial Velocities of F t 

The linearized partial velocities of F i are as follows: 



Vf' =Q , 

(53) 


V 2 =_ ( / l +/ k3)^3 ’ 

(54) 



(55) 

S: F: ( 

1-4 =( 

A +(^3 “4)^2 + [^5 + ^(<72 + Qe)- l W5 -5 ’ 

(56) 

S; F i 

1-5 = 

-l\ [cfc + q l (, ) + / 3 C /6 « 1 + 4 (tf2 + 46 ) ^2 + (4 “ 4? 3 ) £3 > 

(57) 

and 

Vs =4«ia2-('i+4d-'i)a3 • 

(58) 

13. 

. LINEARIZED PARTIAL ANGULAR VELOCITIES 



FOR THE RIGID BODIES OF S 


The following are the linearized partial angular velocities for the system: 


For the zth push-rod: 

S \co P r i = 0 (for r = 1, 2, 3) , 

(59) 


S; P; 

*/®4 = -3 • 

(60) 


S; P; 

1—5 = -2 • 

(61) 

and 

iQ± 6 -q x , 

(62) 

for the zth actuator arm: 
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(63) 


and 


for the flotor: 


and 


S ‘,0>t' = 0 , 


j ; ri: a - / 

/®2* =P\ 


Si A' — f l 
1—3 ~ -3 ’ 


Si m A i - 
1—4 — -3 ’ 


S; , . A; 

/— 5 = —2 


iS; -A; a i 
1—6 ~1\ ’ 


, (1 F _ _„1 
/^1 “ —3 ’ 


Si P A , 1 
/®2 =£i 


5i F -'l 

/^3 -IS 


^1 m F _ pi 
1—4 ~ -3 


_ pi 
/— 5 — -2 


Si p aJ 

/— 6 — 2i ’ 


=0 (r = l 6; i = 2, ..., 8) 


(64) 

(65) 

(66) 

(67) 

(68) 

(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 
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14. LINEARIZED ANGULAR ACCELERATIONS FOR THE RIGID BODIES OF S 

14.1 Linearized Angular Acceleration of Actuator Push-Rod 

S )a P ‘ =u 1 (' I p\ + , 4/?2 + • (76) 

14.2 Linearized Angular Acceleration of Actuator Arm A; 

+ Mg + U 5 Q /2 4" ^4 • (77) 

14.3 Linearized Angular Acceleration of the Flotor F 

S \ QL F = (“2 + + ^ 5/2 + ( _ “1 + “3 + 1{ 4 )/ 3 • ( 78 ) 

15. CONTRIBUTIONS TO THE SET OF GENERALIZED ACTIVE FORCES 
DUE TO THE RIGID BODIES OF S 

15.1 Contributions Due to the Flotor 

On orbit (i.e., neglecting the effects of gravity), the flotor is acted upon by forces and moments 
due to each Lorentz coil, actuator arm, and umbilical, and by direct disturbances. 

c* c* 

Let -F l and -M > represent, respectively, the force and moment exerted by the i th Lorentz coil 
(located at A 3 ) on the flotor, 

where 

F c i=F Ci a i, (79) 

assumed to act at point F ( -, and 

M c i = r F i A i x F Ci = -F C ‘ (l{ + 1\)^ . (80) 
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Let F F ‘ and M F ‘ represent, respectively, the force and moment exerted by the zth actuator arm on 
the flotor, at the zth cross-flexure. Since F F ‘ is a noncontributing force, it can be ignored in the analysis. 
The total moment M F due to the eight cross-flexure springs has value 



(81) 


where k\ is the zth cross-flexure spring stiffness. 

Let F u and M u represent, respectively, the force and moment applied to the flotor by the umbilical, 
where the force is assumed to act at flotor-fixed point F u . Umbilical force F u is given by the equation 

A. A. 

£ = (—*!*! - c x x i)| } + {-k 2 x 2 - c 2 x 2 )§_ 2 

+ (-k 3 x 3 -c 3 x 3 )i 3 +F b ; (82) 

A, 

where q. is some appropriate stator-fixed coordinate system; Xi, x 2 , and x 3 are the umbilical elongations 
in the respective qXi = 1, 2, 3) directions; F_ b is the umbilical bias force in the home position; k l , k 2 , 
and k 3 are umbilical spring stiffnesses; and c x , c 2 , and c 3 are umbilical damping constants. Umbilical 
moment M u is given by 

M U ={-K X 0i-n0\)l x +{-K 2 (p 2 -j 2 ^ 2 )l 2 +{-K 3 (p 3 -Y 3 ^ 3 )l ?) +M b , (83) 


where (j) x , (p 2 , and (p 3 are components of the umbilical angle of twist (j) in the respective q.{i = 1, 2, 3) 
directions; M b is the umbilical bias moment in the home position; K 1? k 2 , and k 3 are torsional umbilical 
spring stiffnesses; and y l5 y 2 , and y 3 are torsional umbilical damping constants. 

Let F°and M /; represent, respectively, the unknown disturbance force and moment acting on 
the flotor. Assume F D to act through the flotor mass center F*. Define/^ and M F to be the zth compo- 
nents, respectively, of F D and M D , componentiated in F x . 

In terms of the above, the flotor’s contribution to the set of generalized active forces, for the zth 
generalized speed, is 
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{Qr= S )v r F “ -F U + S )v r F -F D + S )(0 F r ■ M U +M°+ X k[q[f\ 

V 1=1 


- X ( S )Yr Fi -F Ci + S )(0 F r -M Ci 
i = 1 ' 


The umbilical contributions to the { Q F ’s, namely, S )v r F « ■ F u + S )co F ■ M u , are addressed in the following 

*7 

two sections. The remaining terms of the t Q r ’s are as follows: 

S] l Yf l D =f 2F l D ~f\ [ F 2 D , ( 85 ) 

S )(o{ ■ M D = -M$ , (86) 


s U?{-fC‘) = 0, 

s \o£ \-M c ') = -F Ci {l[ + l\) , 


Si„F ( I uij P ) _ ulj u2J2 , ,3 3 , ,4 4 ,5 5 


/ <T?1 • X hq x f --k\ c l\-k\q\ + k\ C]\ + k\ q\ -k\q\ , 
v'=l ~ 5 ) 


s )y.2 F 1 ' - .'jf./i' -Af'F + l7 F( v hl + v 2 -/iVi) 

S )w^M D = MP+(ql-q\)MP , 

S )t2 (~E- C ‘) = 0 . 


s )e>2 (-M Ci ) = o , 


Vb Ztbi/j =o . 

\i = 1 / 

Vf ' = H + lh\ )Fi D + (v{ + l\q\)F< f , 


S )o)%-M d =M f , 
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(97) 


s )4‘ ■(-E Ci ) = -f C ‘ l l . 

S }sf(-M C ') = F C ‘(li+li) . 

/ g . . ^ \ 

S ) 01 F ■ X </{ / 3 = - &i ' _ *? </f + k\ ql + ki cj\ - k\ qf , 

vi=l > 

, S | / 1-4 • E ° = -4 _l/ 2 + ^(<72 + % ) + l \ q \ F \° 

+ /' (c/3 - q\ j + v} + l\q\ + V3 q\ F® 

+ -v\[q\ + ^6 ) + ( v 2 - 4)<?5 F 3° ’ 

S \(0 F A ■ K D = (q\ + q\)M% - q\M? + M® , 

^4 •(-£ Cl ) = - FQ ( / l- / 3) ’ 

Vf ■(— M c ') = Vf •(-M c ‘) , 

/ « ■ \ ( Q ■ \ 

Si F v- 1 1 l k l Si f , ill 1 
)m ■ X *ltfl/ 3 = /®3 ' X ^1^1/3 , 
v/=l ~ / V/=l / 

5 /vf ■E D = 4^6 + v 2 (<?2 + <76 ) + v 3 ^1° 

+ _v l(<?2 + < ?6) + v 3(ft 1 ~q\) F 2 
+ -v} + v\q\ - f\q\ F ® , 

S \C0 F • M d = -[q\ - q^)M? + M f - + q^M? , 

■ (-M C/ ) = -F c ‘ (/{ + + 4) . 


(98) 

(99) 

( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 

(107) 

(108) 
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( 109 ) 




‘ S| /li6 F°- ^ (91 “ 4;')^| D + [im - v \) F 2 + (4 + v \ c l\ + v 2 - f\ c l\ ) F f > (HO) 

S )0±6 ■ M D = S )G)% ■ M d , (111) 

S Jv?-(-£ Ci ) = ° > ( 112 ) 

S \0±6 \-M c i) = 0 , (113) 

and 

/ « . . • / e . . . \ 

5i ,.F \r jl l Si ,.F v 1 7 / / £ l 

W 6- X*ltfl/ 3 = /®2- I*l«l/ 3 • (114) 

V/=l ~ 7 V/=l / 


Notice the coupling between the unknown-disturbance measure numbers and the generalized coordinates. 
This coupling will make the disturbance input matrix E (in eq. (203)) time-varying. 

15.2 Umbilical Force F u 

Equation (82) expresses umbilical force F u in terms of umbilical-elongation components x 1? x 2 , 
and x 3 , and their time derivatives. These items must be reexpressed in terms of the generalized coordi- 
nates and generalized speeds. 

If the umbilical attachment point F u is at F uh in the home position, then 

Xj = ( r S “ F “ - r s u F uh yi_., ( i = 1, 2, 3) . (1 15) 

But 

/-S(du — I'^u^'uh = fl'idJ'u — fS\F\ _|_ _ j-S\S u _ fS u F u fo , (116) 

where the right-hand- side terms can be expressed by 

r s ^ =llp\-l\a\-l}a\ , (117) 
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(118) 


r_ F i Fu = x Fu f\ + y Fu f l 2 + z Fu f 3 


C C A 1 A, 1 A 1 

r 1 “ = x SuM + ySu s -2 +z Sul3 ’ 
and 

r SuF Kh =x 0 |J + v 0 ^+z 0 ^ , 


for appropriately defined coefficients. 

Now, define the following rotation matrix: 


fl 


'll 

02 

03 

?! 

^1 

I 2 

> = 

r 2l 

r 22 

r 23 

^1 



31 

r 32 

r 33_ 

[-3J 


In terms of the q. coordinate system, equation (116) can now be written in linearized form as 

rV .^, 1 ,+ 4+4 , 

where 

-q\ + <24) • 

- 1} ){<12 + % 

for 

Q - X F U - X S U ~ l 2~ x 0 ’ 

c 2 = x f u - ys u +l l ~ l \ ~ yo > 

and 

C 3 =Z F U ~ Z S U ~ z 0 ■ 



*1 


01 

02 

03 

x 2 

> — 

r 21 

r 22 

r 23 

x 3_ 


/31 

'32 

r 33_ 


C 1 + yF u (<?1 ~ ( 13~ < 14) + Z l 
C 2 ~ Z F U (<?2 + Ql) ~ X F U (dl 
c 3+( l 2 ~ x F u )q\+{yF u 


(119) 

(120) 

(121) 

( 122 ) 

(123) 

(124) 

(125) 

(126) 
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Differentiating, 


Xi 


01 

r l2 

r 13 

x 2 

> — 

r 2l 

r 22 

r 23 

x 3 _ 


_ r 3l 

r 32 

r 33_ 


yp u («1 - «3 - 4 ) + Z F U 4 - l b\ + A 1 («3 + 4 ) 
~ Z F U (4 + 4 ) _ X F U (4 - 4 - 4 ) - 4 (4 + 4 ) 
(^2 _ X F U )4 + (.V/' H - 1 { )(4 + 4) + 44 


(127) 


Using equations (82) and (123)-(127), a linearized expression could now be written straightforwardly for 
the umbilical force F u . 


15.3 Umbilical Moment Mb 


Equation (83) expresses umbilical moment M u in tenus of angle-of-twist components 0 |5 0 2 , 
and 0 3 , and their time derivatives. These items must be reexpressed in tenus of the generalized coordi- 
nates and generalized speeds. 

Let 0 represent the rotation of the flotor, relative to the stator, from the home position, is 
the rotation axis, and 0 is the angle of twist about that axis. Note that 


0i =0«0-| / O‘ = 1 » 2,3) . 


Express h ^ as 


Define rotation matrix Q by 


tin 


~1 ri 

■gif i+g2f 2+83 f 3 



(128) 


(129) 


(130) 
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The linearized 3x3 rotation matrix t Q has elements / Q (/ defined as follows: 


/ 1 
L\ 

- 

'Li 

4-3 



1 

1 1 1 
4 i ~ 43 “ 44 
1 

45 


' 4 l + 43 + 44 


1 

-42 


1 

■46 


“45 
1 , 1 
42 +46 

1 



^1 


*1 


^1 

< 

-2 




[-3 J 


For small 0, it can be shown 13 that 


0 £3 

~83 0 

52 ~5l 


-52 

51 

0 




(l+tr jQ) 


1/2 


(131) 


(132) 


where the post-superscript T indicates matrix transposition and tr jQ represents the trace of { Q. Substitu- 
tion from equation (131) into equation (132), and simplification, yields 


and 


8l=j{q\+q\) , 

1 i 

52 — 45 • 

53 = "^ '(4l - 43“44) • 


(133) 

(134) 

(135) 


Substituting from equations (133)— (135) into equation (129), and transforming into the s } coordinate sys- 
tem by use of ; Q , one obtains the following expression for the spin axis: 


- 1 17 1 i\-i 1 - 1(1 1 l\-i 

»4> = ^42 + 46 J^l + 45^2 “ (4l " 43 " 44 J^3 


(136) 


Since [ n ^ has unit length, 


/0 = 


(42 +46) + (45) +(41 


43 



(137) 
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Use of equations (128) and (136) leads to the following linearized forms for angular position and rotation 
rate: 








/ 1 . 1\ 

9l 


'll 

r n 

63 


^2+%J 

(f>2 

► = 

r 2\ 

r 22 

r 23 

< 

/ * „ \ 

9 3. 


/31 

r 32 

r 33_ 


1 

457 

1 

>«! 

1 

►Si 
4^' — 








/ 1 . 1\ 

4»i 


'll 

r \2 

T3 


\ u 2 + u 6 j 

«5 

92 

► = 

r 21 

'22 

'23 

< 

93. 


731 

'32 

r 33_ 


1 1 1 1 

-1 Hj - M3 - u 4 I 


(138) 


(139) 


From equations (83), (138), and (139), a linearized expression could now be written straightfor- 
wardly for the umbilical moment M u . The flotor’s contribution to the set of generalized active forces 
for the rth generalized speed could then be found by substituting the expressions for F u (sec. 15.2) and 
M u into equation (84). 


15.4 Contributions Due to the Actuator Arms 

The forces and moments acting on the ith actuator arm are due to the respective Lorentz coil 
(located at A 3 ), the flotor (through the ith cross flexure), and the respective push-rod (through the upper 
stinger). The coil force F C ‘ is the only contributing force. The contributing loads, in the above indicated 
order, are as follows: 


C C 

Fyi = f ‘ a\ , assumed to act at point F ; , 

(140) 

M c i = r F i A i x F c i = +l l 4 )qi, x F C ‘ aj = -F C ‘ (/{ +l l 4 ] jg i 3 , 

(141) 

-M p i = -%{/' , 

(142) 

M A i =-k l 2 q2p x -%3& ’ 

(143) 


where lh and l\ are pertinent geometric lengths, and k 2 and £3 are pertinent upper-stinger spring stiffnesses. 
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In terms of the above, the contribution for the /th actuator arm to the set of generalized active 
forces for the /th generalized speed is 


t Qr l = S lYp ■ F C ‘ + S [ 0)f \M Al +M C i ~M F i) . 

(144) 

The individual teims of the ’s are as follows: 


S tYi E Ci =0, 

(145) 

s t (of ■ ( M A i + M c i - M F ‘ ) = 0 , 

(146) 

F c i =0 , 

(147) 

S t 0l2 ■ ( M A ‘ + K c i -M F i) = -tiq l 2 , 

(148) 

S lV F i • F c i = F C ‘ l[ , 

(149) 

• ( M A ‘ + M C i - M F ‘ ) = k[q[ - kiq l 3 - F C ‘ (/{ + 4) , 

(150) 

SvJ.£ c ‘=F c, (/[-4) • 

(151) 

■ ( M A ‘ + M C ‘ - M f > ) = -k[q[ - k& 3 - F C ‘ (l\ + 4) , 

(152) 


(153) 

■ ( M A ‘ + M C ‘ -M F i) = F Ci (l[ + l\jqh + 4) , 

(154) 

^■^=0 , 
and 

(155) 

* ' ( M A/ + M c i -M f i) = -k& 2 . 

(156) 


Notice the coupling between the control inputs and the generalized coordinates. This coupling will make 
the disturbance input matrix B (in eq. (204)) time varying. 
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15.5 Contributions Due to the Push-Rods 


The contributing loads on each push-rod are moments and -M A ‘ , where (using pertinent 
lower-stinger stiffnesses) 


M P ‘ = ~k\q\l 3 ~ k W ^-2 ~ k (/l6q\ . 

(157) 

The contribution for the /th push-rod to the set of generalized active forces for the rth generalized 
speed is 

lQp= S l(0 Pi ■ ( M p i - M a i ) . 

(158) 

p 

The individual terms of the l Q r l ’s are as follows: 


S [(Or‘ ■ (M P ‘ - M A ‘ ) = 0 (r = 1, 2, 3) , 

(159) 

^(0% ■ [M p i - M A ‘ ) = k^q‘ 3 - k l 4 c/ 4 , 

(160) 

S ,^-(M P . -M A ‘ ) = -k l 5 q‘ 5 , 
and 

(161) 


(162) 


16. CONTRIBUTIONS TO THE SET OF GENERALIZED INERTIA 
FORCES DUE TO THE RIGID BODIES OF S 


A- / A- 

Represent (by I -1 1 ) the central moment/product of inertia of the /th actuator arm for the j and 

^ * 

k body-fixed coordinate directions q l j and q^. Define push-rod inertias lf^ P ‘ analogously, where the 

F/F* 

single subscript indicates that the axes are assumed to be principal axes. Let I ^ represent the central 


*1 ^1 

inertia scalar of the flotor for the flotor-fixed coordinate directions / . and /. . Use the symbol H to 

— J — * 


represent an angular momentum vector. Associated post- superscripts on H have the same meanings as 
for the inertias. The contributions to the generalized inertia forces for S can now be expressed. 
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16.1 Contributions Due to the Push-Rods 


/ *\ P i 

The contributions I / Q r I to the generalized inertia forces due to the zth push-rod are as follows: 


(,&r 


.S', P. 
■■ \v r ' 


S Ln P i . f P ‘ 


P \ 


m P . \+ iQl r l ■ ~iK 1 =Q (r = 1,2,3; / = 1, ...,8) , (163) 


{iQl) P> = s \ u‘ ■ {- m p t s j a Pi )+1 J • -/ ^ Pi 


mp^Pi) +I ^ P ‘ 




(164) 


(/Gsf = S ‘i4‘ {■ -m Pi Si ia Pi )+'/«? i-iiL 


p. / p* \ p/p* i 

f l/ l = —I pPl W5 


(165) 


(/Cfef = s ',4 (~"‘p, s i^*)+i®6 (-,s. p ‘ ,p ‘ 

>"p {p'-i r +i/ pi 


»6 


(166) 


16.2 Contributions Due to the Actuator Arms 


A: 


The contribution ([Q* j ‘to the generalized inertia forces due to the ith actuator arm for the rth 
generalized speed is 


(iQr) A ‘ = S 0 ■ ( -m A . S \a Ai 


-,H 


A:/A 


*\ 


(167) 


Then the individual, nonzero terms of the 


K) A ' 


s are as follows: 


A A i ( A A* ) \( i \ 2 ■/, 

W ■ | -m A . /« ' ) = - rn A . U a 2 ) u 2 + 


a l[ a 2 +4) 


• l l l -l 
lift — 


(168) 
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16.3 Contributions Due to the Flotor 


The contribution 



to the generalized inertia forces for the rth generalized speed is 




S. f 

l n 1 


a 



( 169 ) 

( 170 ) 

( 171 ) 

( 172 ) 

( 173 ) 

( 174 ) 

( 175 ) 

( 176 ) 

( 177 ) 

( 178 ) 
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/ *\F 

Then the individual terms of the I jQ r ’s are as follows: 


Si F* ( S: f 
)v{ \-m F \a r J = -m F 



4 + f\ 44 


+ 


-fl v l ~ fl v 2 ] 4 ~\fl (4 +l \) + f\ v \ 4 + .444 + /lM4 } > ( 179 ) 


M- -is. 


-,H F / F 


* \ 


t F/F j .1 . / \ T'//-' ./ ,7'//'’ / ./ . .A 

— 7^ | I «2 + ^6 I + 7^2 + 7^ I M3 + M4 — Mj I , 


F/F t F IF 


(180) 


/4 ' ) = -m F {v\f\u\+ (4) + (4 


\2 / 1 \2 


1 1 1 
“2 _V 1 V 3 M 3 — Vi V3W4 


■ 444 + 


(vi) 2 + vi(vi +/j) 


4 


(181) 


M • (-/£ F/F *) = -4T* (4 + 4)- /f 2 /F * 4 - 4T* (4 + 4 - 4) . asg 


S U,F 


)Y 3 \ ->n F S \a F ) = -m F {[-/ 2 ‘ 4 - f\ 4 ] 4 - 4 44 + (4) +(4) 

4) 2 +4(4+^!) 


“3 


+ 




(183) 


S H 


o) F • = -i F { F * (4 + 4)- 4 /F * 4 - 4 /F * (4 + 4 - 4) , a 84) 


S )Y4 ■{- fr >F S 'i c / j = -m/ 7 j-/ 2 l ( / 3 + 4)4 


All „1„M , 
- /l ^1 ll \ - v 1 v 3 m 2 + 


(4f + 4(4+4) 


• 1 

w 3 


+ 


4) 2 + (4 + ^) 2 4 - 4(4 + 4)4 - 444} . (i85) 
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(186) 





~ m F { /2 v 3“l " v l v 2«2 

- V 2 V 3 u\ - V^vJ) + 4 )^ 4 + [( v l)" + ( v 3 ) 


2 




and 



S m F 

1—6 


(-,a F/f *) = -I[r (4 + 4)- I?r 4 - /, f 3 ,F * (4 + 4 - 4) . 

17. EQUATIONS OF MOTION FOR THE SYSTEM 


(187) 

(188) 


(189) 


(190) 


17.1 Kinematical Equations 

There are 48 kinematical equations for the system, one for each generalized speed: 

u) = q), 0 =1,..., 6; / = 1, ..., 8) . (191) 

17.2 Dynamical Equations 

Six dynamical equations are obtained using the following process. First, add the respective 
contributions of the 17 rigid bodies to the set of holonomic generalized active and holonomic general- 
ized inertia forces, for each generalized speed (i.e., r = 1, ..., 48). The holonomic generalized active 
force for the rth generalized speed is 
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(192) 


Fr= l{F r ) j , 

7=1 

where {F r ) 1 is the contribution to the set of holonomic generalized active forces due to the 7 th rigid 
body. That is, 


F 8 A 8 p 

F r = iQr + 1 lQ? + I iQr ■ (193) 

i=l i= 1 

Likewise, the contribution to the set of holonomic generalized inertia forces is 


7=1 




(194) 


Second, develop the relationship between the dependent and the independent generalized speeds in the 
form 


6 

= I 
5=1 


A rs u s +B r (/ = 2, ..., 8 ; j = 1, ..., 6 ; r = 7, ..., 48) 


(195) 


where the factors are the six independent generalized speeds, and the terms u‘j are the 42 dependent 
generalized speeds. A rs and B r are scalars, derived from the nonholonomic constraint equations (see 
sec. 17.3). The nonholonomic and holonomic generalized active forces are related to each other as 
follows: 


48 

F r =F r + I F s A sr (r = 1, ..., 6 ) . 
s=l 


(196) 


Similarly, the nonholonomic and holonomic generalized inertial forces are related to each other as follows: 


48 

F* = F* + X FsAsr (r = 1, .... 6 ) 
s=7 


(197) 


Kane’s dynamical equations,* then, are 


F r +F r = 0 (r = 1 , ..., 6 ) . 


(198) 


* This step was erroneously omitted in references 10 and 11, resulting in incorrect incorporation of system constraints. The authors 
discovered and corrected the error during the course of model verification. 


31 



17.3 Constraint Equations 


The kinematical and dynamical equations together are 54 in number: 48 kinematical and 6 
dynamical. Since the complete set of equations for an //-degree-of- freedom system numbers 2 n, and 
since the system S has 48 degrees of freedom, 42 more equations are needed to describe completely the 
motion of the system. These missing equations are the holonomic constraint equations (in nonholonomic 
form) for the dependent generalized speeds «'• ( i = 2 , ..., 8 ; j = 1 , ... 6 ). 

Since the velocity and angular velocity of the flotor center of mass F* is the same irrespective of 
the actuator path chosen for describing its position, a set of constraint equations can be written in vector 
form using the following: 


and 


Si f* 

l v - 


s >dA F 

dt 


* \ 


\v F 0 = 1, y = 2 ,..., 8 ) 


(199) 


S \Q) F * = Sj l co F \i = \, j = 2 ,. .., 8 ) . 


(200) 


A, 

If one expands equation (199) and resolves them into a common coordinate system (here, the f. coordi- 
nate system), one obtains the following 21 (motion) constraint equations: 


1 

i 

i 

T 

i 

U\ 


Aflj-flfy 

T 

u\ 

-^ 3/2 j + v 2fij 


u 2 


+ v \f3j 


u 2 

-Mi + AAj 


u 3 


-v\f\j + Mj 


u 3 

M-ibAj+AAj 


i 

U4 


~(v\ ~l\)f\j + v\f 2 j 


ll\ 

. 

TT 

l 

> 


u 5 


Mj -Mj 


4 

u l_ 

-V3f2j+(V2+l3)f3j_ 


u l 6 


- v 3f2j+( v 2 + l bf3j 



0 = 2, — , 8 ; 7 = 1, 2, 3) , (201) 


Similarly, if one expands equation (200) and resolves them into a common coordinate system (here again, 
the /. coordinate system), one obtains the remaining 21 (motion) constraint equations: 
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u 3 

V 
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1 

«4 
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m 5 
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_ f) _ 


1 

l u 6\ 


(202) 


(It should be noted here that, although equations (201) and (202) are in nonholonomic form, the constraints 
they represent are actually geometric.) 

18. STATE-SPACE FORM OF THE EQUATIONS OF MOTION 


The time derivatives of the constraint equations can be used with the constraint equations them- 
selves, to write the kinematical and dynamical equations (eqs. (191) and (198), respectively) in the 
following descriptor form: 


/ 0 

m. 

0 

N~ 

j 2 u 

~0~ 

(0+ 

O' 

O M 

b r \ 

K 

c 

[w 7 ] 

B 

E 


{d}. 


(203) 


The state vector consists of the 48 coordinates q and the six independent generalized speeds u 1 , where 


u 


/ = 


1 

U\ 


1 

u 2 


1 

« 3 


1 

«4 


1 

«5 



(204) 


The constant matrices M, K, and C are system mass, stiffness, and damping matrices, respectively. The 
symbols I and O represent, respectively, an identity matrix and a zero matrix of appropriate dimensions; 
vector i contains the eight control currents to the Lorentz coils; and vector d is the disturbance vector. 
The input matrices B and E are time-varying matrix functions of the coordinates. A is a constant matrix 
that incorporates the kinematical equations; N and C together incorporate the holonomic constraints. 

The disturbance term \E\ [d] accounts for the umbilical bias force F_ b and moment M b , and the 
unknown direct disturbance force F D and moment M D . Recall that in the development of the foregoing 
equations the angular acceleration of the stator was assumed to be negligible. However, the translational 
acceleration of the stator, although presumably unknown, cannot be neglected. In fact, that acceleration 
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is the source of the umbilical contribution to flotor g-jitter. To include this indirect disturbance 
contribution, one must add an (unknown) indirect acceleration disturbance term a 1 to each of 
equations (26)-(28). Along with the other disturbances, this indirect disturbance will appeal' in the final 
term of equation (203). 


19. MODEL VERIFICATION 


AUTOLEV™ software, marketed by Online Dynamics, Inc., was used to create a full nonlinear 
model of ARIS including the actuator (rigid-body) dynamics. 14 AUTOLEV was then used to develop 
and verify the linearized equations presented in this paper. 

An independent model was developed using the DENEB Envision software, with current CAD 
models of an ARIS -outfitted ISPR. 14 This model was used as an independent (static) check of the 
actuator kinematics. 

For dynamic validation, various loads were applied to the flotor, and the nonlinear AUTOLEV 
model was allowed to move in simulated response. Motions were permitted which greatly exceeded the 
rattle space constraints (with angles permitted up to -90°), to test thoroughly the nonlinear model. The 
six angles for a single actuator (AUTOLEV model) were inserted into the Envision model, and the 
remaining 42 actuator angles were compared between the two nonlinear models. The respective angles 
were in consistent agreement, even with these large motions, to within less than half a degree. 


20. FUTURE WORK 


The next tasks will be the addition of umbilical forces to the AUTOLEV equations, the imple- 
mentation of a linearized model in MATLAB®, and MATLAB model verification. The verification 
procedure will first entail comparing the eight position vectors from a common point on the stator to the 
flotor center of mass, as traced through the eight actuators, with the flotor centered in its home position. 
Then the procedure will be repeated with the flotor moved statically from its home position, in six 
degrees of freedom. 

For dynamic verification, various loads will be applied to the flotor to verify that the eight 
position vectors track for motion inside the rattlespace (i.e., small angles). Simulations of the nonlinear 
AUTOLEV model will be compared with simulations of linearized system models (one with and one 
without actuator dynamics) to determine the simplest model suitable for controller design. 

System dynamics will be incorporated into the Envision model, along with the capability of 
state-space, discrete time control. The MATLAB and Envision models will then be available, respectively, 
for centralized, state-space/optimal controller design and for closed-loop system simulation. 
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